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Abstract 

We investigate the influence of the specimen velocity on the magnetic flux 
leakage with the aim of selecting the optimum sensor locations. Parametric 
numerical simulations where the specimen velocity was in the range [0.1-20] 
m-s“^ were carried out. As the specimen velocity is increased, the magnetic 
held varies from being symmetrical to being asymmetric. For the radial mag¬ 
netic induction, the peak to peak value moves from the centre of the bridge 
towards the direction of the specimen movement. For the axial magnetic in¬ 
duction, the specimen velocity influence is dependent on the sensor location 
and a signal-velocity independent region was discussed. 

Keywords: Nondestructive testing, magnetic flux leakage, specimen 
velocity, peak to peak value. 


1. Introduction 

Effective nondestructive testing (NDT) can prevent disasters similar to 
the Buncefleld incident jl|, resulting in an estimated loss of £894 million. A 
great number of NDT methods have been developed. These include: visual, 
radiographic, ultrasonic, eddy current, thermal infrared, acoustic emission, 
and magnetic particle testing, etc.. In this paper, we focus on the magnetic 
flux leakage (MFL) method, which originates from the magnetic particle 
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technique. The principle of the MFL method can be understood as follows: 
when a magnetic field is applied to a ferromagnetic material, the leakage of 
the magnetic field which is caused by the discontinuity of the geometry can 
be captured by the sensors, such as Hall probe etc.. The main reason for the 
leakage is due to the difference, ~ C>(10), of the magnetic permeability of the 
mediums at the interface j^. The first application of MFL can be dated back 
to 1868, the British navy used a compass to examine a magnetised cannon 
to search for defects j^. This technique was extensively applied due to the 
development of magnetisation techniques for examining defects in pipelines, 
pressure vessels, and wheels, etc. in the 1960s. The defect characteristics, 
such as shapes, dimensions and locations, can be determined by the leakage 
signals and a large amount of relevant numerical and experimental research 
has been carried out. Practically, the nature of such NDT problems is tran¬ 
sient rather than static and simply using the results obtained from static 
simulations to predict transient problems may cause errors, especially under 
the condition of high velocity. A velocity induced current can be generated 
by a conducting material moving in a magnetic field jj and this phenomenon 
can alter the distribution of the magnetic field jEl, @,13, 0] . The distortion of 
the magnetic field points in the specimen movement direction and is indepen¬ 
dent to the orientation of the magnetizing source. This signal deformation 
can influence the efficiency of the NDT, especially for determining the defect 
location. With the aim of compensating for this leakage signal deformation, a 
scheme was validated against experimental results 0]. However, the optimum 
location of the sensor is still not well defined, especially for high velocities. 
In order to cover an extensive range of parameter in term of specimen ve¬ 
locity, we shall combine a detailed analysis of a pure two-dimensional (2D) 
geometry with tarceted numerical simulations of the full three-dimensional 
(3D) problem BSQ- Clearly, this is the task which shall be undertaken 
in this paper. 

We shall first describe the configuration and numerical setup in section [21 
In section |21 we investigate the distortion of the magnetic field due to the 
specimen velocity effect. We then discuss the influence of the specimen ve¬ 
locity on the radial magnetic induction By in section 13.11 and axial magnetic 
induction in section 13.21 General conclusions are summarized in section 
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2. Configuration and numerical setup 

2.1. Configuration 

We consider a ferromagnetic specimen (conductivity a = 6.993 x 10® 
S-m“^, having B-H curve) moving with velocity (in m-s“^ in this paper) 
under a magnetic flux leakage evaluation system which uses permanent mag¬ 
nets as the magnetising source as shown in FigO The direction of specimen 
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Figure 1: Schematic representation of numerical simulation model for MFL evaluation 
system. All dimensions are in mm. The materials of bridge and brush are identical to the 
specimen. Dru is defined as the relative distance between and the centre of the defect. 
The left edge (resp. right edge) of the defect is defined as front edge (resp. back edge). 

movement is taken as the x-axis. A near-side rectangular defect is located on 
the specimen and the reason for the selection of a rectangular shaped defect 
is because 2D finite element methods (FEM) can movide sufficient informa¬ 
tion for the sharp-shaped defect characterisation [uj . The dimensions of the 
bridge, magnet, brush and the defect are shown in the figure as well. 

Sj denotes the potential sensor locations where i G [1,11] (along +x axis) 
and j G [1,3] (along +y axis), where i and j are the column and row num¬ 
bers. These potential sensor locations Sj are equally spaced 3 mm apart. 
Furthermore, the relative distance between Sj and the centre of the defect is 
defined as Dri. Dri > 0 {resp. Dri < 0) denotes the centre of the defect is 
approaching {resp. is departing away from) the sensor. 
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2.2. Governing equations and numerieal setup 
2.2.1. Governing equations 

We study the 2D magnetic flux leakage problem using Ansoft Maxwell 
(version 14.0) FEM software. The governing: equations for the 2D transient 
MFL problem can be expressed as |l2l |: 

1 

V X —V X A = Js — cr—- aW + V x He + aVg x V x A, (1) 

p at 

B = V X A, (2) 

where /U, A, J^, a, V, and He are the permeability, magnetic vector 
potential, sonree current density, electric conductivity, electric potential, ve¬ 
locity and the coercivity of the permanent magnets, respectively. For the 
MFL evaluation system investigated in this paper, Equation [T] becomes: 


V X 


1 

— V X A = —a ——h V X He + erVs x V x A, 
fi at 


( 3 ) 


which is solved using FEM with infinite boundary conditions, together with 
Eqnation [21 The whole computing domain is discretised into 2D triangular 
elements, and the first-order implicit Euler method is adopted as the time 
discretization method. A computational domain percentage, 200% US , is 
selected. 


2.2.2. Numerical setup 

We first performed a mesh sensitivity analysis to determine the mesh 
density required for the simulations to give a mesh independent solution. The 
main characteristics of the meshes tested are provided in Table [H By and 
denote the radial and axial magnetic induction. We compute the errors cbx 
(respectively esy) on B^ (respectively By) relative to the finest mesh M3 at 
different sensor locations, and denote the discrepancies of B^ and By 

between different meshes at the sensor point Sj: e.g. = |1~ 1; where 

k=l and 2. Both these decrease when the total number of mesh elements 
increases, which shows good convergence. Therefore in this paper, all the 
simulations are based on mesh M2, which ensures a satisfactory precision at 
a reasonable computational cost. For M2, the length of element size 0.5 mm 
for the bridge, magnet, brush and the specimen. 

A further validation test was carried out to determine the required setup for 
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Table 1: Main characteristics of the different meshes and errors in the magnetic induction 
for Vs=5 m/s. 


Meshes 

Ml 

M2 

M3 

Number of nodes between two brushes 

2526 

10322 

39954 

Total number of nodes 

2.1x10^ 

8.1 xlO^ 

1.1x103 

4 

1.3x10-3 

3.5x10-^ 

- 

4 

2.2x10-2 

3.1x10-3 

- 

4 

9.6x10-^ 

3.2x10-^ 

- 

^By 

1.3x10-2 

3.0x10-3 

- 


1.2x10-3 

3.6x10-^ 

- 

By 

2.3x10-2 

3.2x10-3 

- 


the transient simulations. This was done by comparing the results obtained 
from the static solver and the results obtained by using the transient solver 
with no velocity (Vs=0). We compare at different sensor points (Sj, i=l, 
3, 5 and 6) at Dri = 0 and the results show that the discrepancies of are 
4.9x10“^, 3.5x10“^, 2.5x10“^ and 2.1x10“^, respectively. Therefore, the 
numerical setup is valid. 

3. Results 

Compared to the static simulations, one of the most important influences 
of specimen velocity is the distortion of the magnetic field, as shown in Figj2l 
The snapshots are the distribution of the magnetic flux lines and are captured 
at the moment when Dri = 0 for Sg at different specimen velocity (Dri was 
defined in Figil]). For the low values of the specimen velocity, the magnetic 
field is symmetric to the centre of the defect, as shown in FigJ21 (a)- As Vg 
increases (b - d), the distortion of the magnetic field occurs. The distortion 
here introduces the asymmetric feature with respect to Sg. This distortion 
is mainly caused by the velocity induced current, which was well discussed 
013, 13. [f^ . In this paper, we mainly focus on the phenomenon which is 
resulting from this distortion. 

3.1. Radial magnetic induction By 

3.1.1. General feature 

The distribution of By with Dri at Sg for different specimen velocity is 
shown in FigEl The result shows that the minimum value of By (respectively 
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Figure 2 : Snapshots of distributions of magnetic flux lines when Dri = 0 at 85 . From top 
to bottom: = 0.1,1,10 and 20. The maximum value is -0.0144 (in red) and minimum 

value is -0.0146 (in blue) Wb/m. The magnetic field distortion occurs as is increased. 
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Figure 3: Distribution of By vs. Dri at Sg for different Vs; 0.1 (□), 5 (o), 10 (A), 20 
(v)- The maximum and minimum values are presented when the sensor meets the defect 
edges. By moves downwards as Vg is increased. 

maximum values of By) 5™” (respectively 5™“*) occurs near the front (re¬ 
spectively back) edge of the defect for all values of specimen velocity. This 
can be understood as follows: in the vicinity of Dh= 0.005, sensor Sg meets 
the front edge of the defect. The leakage has a trend to return into the spec¬ 
imen due to the high permeability of the specimen. Under this condition, an 
angle between the x-axis and B is presented, which has the effect of enlarging 
the negative component of By. When =-0.005, Sg meets the back edge of 
the defect. The leakage occurs due to the existence of the forthcoming defect 
and this results in By reaching B^^^. Furthermore, distribution of By on Sj, 
where i G [1,11] and j G [1,3], follows the same trend. This indicates that 
Vs will not influence the effect of the defect edges on By, compared to the 
static situation. 

The distribution of By shifts downwards, as specimen velocity is increased. 
As we discussed before, the distortion of the magnetic field occurs at high 
value of the specimen velocity. This distortion causes the anticlockwise ro¬ 
tation of the magnetic induction between the two brushes of the MFL eval¬ 
uation system, as shown in FigJU Take the magnetic induction at Sg as an 
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Figure 4: Distribution of magnetic induction B and iso-lines of By at Dri=0 for S@ for 
Vs=0.1 {top) and 20 {bottom). The magnetic induction B rotates anticlockwise as 
increases. 


example: at low values of Vg, B is almost parallel to the a:-axis for the sensor 
location. The radial component is almost zero. However, at highs value of 
Vs, B rotates anticlockwise and an angle between B and a:-axis occurs. This 
angle results in the non-zero magnitude of By. We hnd that the variation of 
By obtained at Sg is proportional to the specimen velocity (maximum value 
is 20 m-s“^). 


3.1.2. Peak to peak value B^ 

In practice, the sensitivity of the MFL signal is not only dependent on 
it’s magnitude, but also on the variations of the signal jl^. We then focus 
on the variations of the radial leakage signal. The peak to peak value for By 
is dehned as follows: 


_ ^max 


TDrmn 


(4) 


The distribution of B^ at different Si for different is shown in FigJSl 
Firstly, the maximum B^ occurs at Sg when the specimen velocity is low, 
e.g. Vs =0.1 (□ in the figure), which is located centrally between the two 
permanent magnets, for all values of j. For low specimen velocities, the 
magnetic field is symmetric about the centre of the defect at Dj-i = 0 and 
Bmax mainly caused by the edges of the defect. B^ varies 

a little at Sj and the curve remains relative flat. Secondly, as specimen 
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Figure 5: Variations with Sj for different Vg. Left: j=l; middle: j=2; right: j=3 for 
Vs=0.1 (□), 5 (o), 10 (A), 20 (v)- BP moves towards the specimen movement direction 
as Vg is increased. 

velocity is increased, the maximum value of occurs further towards the 
direction of the specimen movement. This phenomenon is mentioned by Shin 
jsl, whose hndings are based on an MFL system which uses direct current 
electromagnets as the magnetising, however, the underlying reason for the 
phenomenon was not well discussed. The phenomenon is mainly due to the 
distortion of the magnetic field caused by specimen movement. This can be 
understood as follows. As we mentioned before, and occur at the 
moment when the sensor meets the front edge and back edge, respectively. 
FigEl shows the By distribution at the moment when the sensor meets the 
front and back edge for the sensor S4, Se and Ss at Vg=10. The results show 
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Figure 6: Snapshots of By distribution when 84 , Se and Sg meet the front edge {top) and 
back edge (bottom). Top: the maximum value is -0.058 (in red) and the minimum value 
is -0.079 (in blue). Bottom: the maximum value is 0.03 (in red) and the minimum value 
is 0.0045 (in blue). 

that from Sg to S 4 , for the magnitudes of By decreases from 0.079 

to 0.058 (26.6%). However, for the magnitudes increases from 0.0045 


9 

















































to 0.03 (566.7%). This unbalanced increase results in moving forwards 
in the specimen movement direction. For high values of Vg, the optimum 
sensor location is away from the middle of the bridge in the direction of V*. 

3.2. Axial magnetic induction 
3.2.1. General features 

The distribution of axial magnetic induction with D^i at the sensor lo¬ 
cation Sg for different specimen velocity is presented in FigjTl The results 



Figure 7: Distribution of Bx vs. Dri at Sg for different V^: 0.1 (□), 5 (o), 10 (A), 20 (v)- 
For Dri=0 mm, Bx decreases as Vs is increased. 

show that for all values of V*, the maximum magnitude of B^ occurs in the 
vicinity of Dri=0. This is due to the fact that at that moment B is almost 
parallel to the x-axis. This helps B^ reach its maximum magnitude, as shown 
in the figure. The variation of B^ at all potential sensor locations Sj follows 
the same trend. As is increased from 0.1 to 20, the maximum magnitude 
decreases by 11.7% in the present research. 

3.2.2. Bx on Sfi for selected Dri 

In the section 13.2.11 we focus on the sensor point Sg. In this section, we 
turn our attention to the B^ distribution on all the sensor locations. 

As the lift-off value (j) is increased, the magnitudes of B^ decrease (FiglH] 
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left). This simply reflects the influence of the sensor lift-off values. In order 
to capture the leakage signal effectively, the sensor should be located near 
the specimen. However, this phenomenon is weaken as the specimen velocity 
is increased (Fig^right). 



Figure 8: Left: at different sensor locations for j=l and 2 at Vg 


10. Right: 




Sir 


at different sensor locations for different specimen velocity at =0.005 and 0. 


The variation of with Si at different specimen velocity for different Dj-i 
(0.005, 0 and -0.005) is shown in FigjOl In the figure, from top to bottom, 
2 and 3, respectively and from left to right, Dri=0.00f>, 0 and -0.005, 
respectively. For Vs=l (□ in the figure), the asymmetric (with respect to 
Bx on Sg) feature occurs for the selected at j=l. The sensor which 
first meets the defect measures the strongest leakage signal of B^. As j is 
increased, this feature becomes insignificant. This reflects that increasing j 
can reduce the influence of Vs, however, the strength of the Bx signal is also 
decreased. As is increased, this asymmetric distribution feature becomes 
more dominant, even for ^’=3. 

Interestingly, for the geometry we investigated the results show that the 
influence of on the magnitude of Bx is highly dependent on the sensor 
location for all values of Dri- At Dri=0.00b and j=l, for the sensor locations 
S} i G[1,5], the magnitude of Bx decreases as V* is increased, and conversely, 
for S} i g[7,11], the magnitude of Bx increases as is decreased. The 
velocity has little influence at Sg under this condition. Here, we define Sg 
as the leakage-velocity independent region. The existence of this region is 
mainly due to the differences between Bx gradients for different V^. This 
location moves in the direction of V* as j is increased, for all selected D^i- 
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Figure 9: Variations of Bx with at Dri=0.005, 0 and -0.005 for different Vg: 1 (□), 5 
(o), 10 (A), and 20 (v)- From top to bottom: j=l, 2 and 3; from left to right: Dri=0.005^ 
0 and -0.005. Leakage-velocity independent region is highlighted with a red square. 

In practice, it is helpful to locate the sensor in the vicinity of this region to 
decrease the influence of from the axial magnetic held point of view. 

4. Conclusions and future work 

We have conducted a detailed analysis of the influence of specimen ve¬ 
locity on the magnetic flux leakage signal using finite element analysis. The 
main results can be summarised as follows: 

• deformation of the magnetic held occurs as the specimen velocity is 
increased, which is due to the velocity induced eddy current effect. 

• For the radial magnetic induction By^ the maximum variation of the 
leakage signal moves from the centre location between the two brushes 
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towards the specimen movement direction. The reason for this move¬ 
ment is that the leakage that escapes from the specimen increases much 
faster than the leakage which returns back to the specimen as the spec¬ 
imen velocity is increased. This movement becomes insignihcant as the 
sensor location is further away from the specimen. This indicates that 
the optimal sensor location is not in the middle of the bridge, but at a 
location some distance from this point in the direction of the specimen 
movement, especially for higher speed MFL evaluations. 

• For the axial magnetic induction Bx, the influence of the specimen 
velocity is highly dependent on the sensor location, for the geometry we 
investigated. A potential sensor location which can capture a constant 
value of Bx leakage at different specimen velocities is fonnd for the 
geometry we investigated. This gives us greater flexibility to select 
suitable sensor locations. 

Future work will focus on 3D simnlations and experimentally investigating 
the optimum sensor location for high speed measurement based on peak to 
peak valued. Furthermore, the leakage-velocity independent region which is 
based on the axial magnetic induction will be investigated experimentally to 
confirm the simnlation results presented in this paper. 
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